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In 1876, Edouard Lucas discovered a method for proving that a number is 
prime or composite without searching its factors. His method was based 
on the properties of the Lucas Sequences. He first used his method for 
Mersenne numbers and proved that 2^^^ — 1 is a prime. In 1930, Derrick 
Lehmer provided a complete and clean proof. This test of primality for 
Mersenne numbers is now known as: Lucas-Lehmer Test (LLT). 

Few people know that Lucas also used his method for proving that a Fermat 
number is prime or composite, still with an unclear proof. He used his 
method for proving that 2^^ + 1 is composite. Lehmer did not provide a 
proof of Lucas' method for Fermat numbers. 

This paper provides a proof of a LLT-like test for Fermat numbers, based 
on the properties of Lucas Sequences and based on the method of Lehmer. 
The seed (the starting value 5*0 of the {Si} sequence) used here is 5, though 
Lucas used 6. 

Primality tests for special numbers are classified into — 1 and + 1 
categories, meaning that the numbers — 1 or + 1 can be completely or 
partially factored. Since many books talk about the LLT only in the N 
chapter for Mersenne numbers A'^ = 2"^ — 1, it seemed useful to remind that 
the LLT can also be used for numbers A^ such that A^ — 1 is easy to factor, 
like Fermat numbers A^ = 2'^" + 1, by providing a proof dla Lehmer. 

Theorem 1 

Fn = 2^" -\- 1 (n ^ 1) is a prime if and only if it divides 5'2"-2 , where 
So = 5 and Si = S^^_-^ - 2 for i = 1, 2, 3, ... 2" - 2 . 

The proof is based on chapters 4 (The Lucas Functions) and 8.4 (The Lehmer 
Functions) of the book "Edouard Lucas and Primality Testing" of H. C. 
Williams (A Wiley- Interscience publication, 1998). 

Chapter 1 explains how the (P, Q) parameters have been found. Then Chap- 
ter 2 provides the Lehmer theorems used for the proof. Then Chapter 3 and 
4 provide the proof for: F„ prime =^ Fn \ 5'2n„2 and the converse, prov- 
ing theorem 1. Chapter 5 provides numerical examples. The appendix in 
Chapter 6 provides first values of Un and Vn plus some properties. 



AMS Classification: 11A51 (Primality), 11B39 (Lucas Sequences), 11-03 
(Historical), 01A55 (19th century), 01A60 (20th century). 



1 Lucas Sequence with P = 

Let So = 5 and Si = Sf^^ - 2 . 5i = 23, ^2 = 527 = 17 x 31, ... 

It has been checked that: | ^'f^~'^ ~ ^ (™od -F„) for n ^•••^ 

52n_2 ^ (mod Fn) for n = 5.. .14 

Here after, wc search a Lucas Sequence {Um) m>o and its companion (V^) m^o 
with (P, Q) that fit with the values of the Si sequence. 

We define the Lucas Sequence Vm such that: 

y-2k+i = Sk (1) 



Thus we have: 



' V2 = So = 5 
V4 = Si = 23 
^8 = ^2 = 527 



V4 = Vi - 2Q2 



Vs = Vi - 2Q 



If (4.2.7) page 74 ( V2n = - 2Q^ ) applies, we have: 

and thus: Q=^/^=y^ = ±l. 

With (4.1.3) page 70 ( T^+i = PVn - QVn-\ ), and with: 

' Vb = 2 
< Vi=P 

V2 = PVi - QVo = P2 _ 2Q 

we have: P = VWT2Q = or . 

In the following we consider: (P, Q) = [s/l, 1) . 

As explained by Williams page 196, "all of the identity relations [Lucas 
functions] given in (4.2) continue to hold, as these are true quite without 
regard as to whether P, Q are integers" . 

So, like Lehmer, we define P = \/P such that R = 7 and Q = 1 are coprime 
integers and we define (Property (8.4.1) page 196): 



Vn when 2 I n _ f Un/VR when 2 I n 

r- Vn = \ 

Vn/yR when2-|'n [ Un when 2 \n 



in such a way that Vn and C/„ are always integers. 

Tables 1 to 5 give values of Ui , Vi , Ui (mod P„) , Vi (mod P„) , with 
(P,g) = {V7,l) , forn = 1,2,3,4 . 
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2 Lehmer theorems 

Like Lehmer, let define the symbols (where (Vb) Legendre symbol): 

' e = e{p) = (D/p) 
< a = a{p) = {%) 

. r = rip) = {%) 
The 2 following formulas (from page 77) will help proving properties: 

(4.2.28) 2"'-'Umn = E (si + lj^'^'''''^""^''^'^ 



i=0 
m I 



(4.2.29) 2— = ^ (^^J^*?72V--2^ 



Property (8.4.2) page 196 : ^ 

(U,^{%) (modp) 
If p is an odd prime and pi Q, then: < 

1f,^(R/p) (modp) 

Proof: 

Since p is a prime, and by Fermat little theorem, we have: 2^"^ = 1 (mod p). 

• By (4.2.28), with m = p and n = 1, since Ui = 1 and V\ = P, we have: 

p-i 

2f-iC/^ = (^^P^-' + (^^DPP-^ + ... + (^^d'^P^ 

Since (^) = (mod p) when < i < p and (p = 1 , we have: 
C/p = I7p = = (D/p) (modp) 

• By (4.2.29), with m = p and n = 1, since ?7i = 1 and Vi = P, we have: 

2-'F, = (;;) P' + Q Di-- + ... + ^ J 

Since (JJ) = 1 , and (^) = (mod p) when < i < p , we have: 
Vp = PP and Vp = PP-^ = = (%) (mod p) 

□ 
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Property (8.4.3) page 197 : 

p odd prime and p \ Q =^ p \ Up-ae 

Proof 

By (4.2.28) with n = 1, Vi = P, since p is a prime and {R, Q) = 1, we have: 

• With: m = p+1 

i=0 ^ 

VUp+i = ip+ l)PP + {p+ l)p[...] +{p+ l)D^P + OD^P-^ 
2PUp+i = PP + D"-2^P + p[...] =P[(P2)^+D^] 

= 2PUp+i = r'^+d'^ = {%) + (D/p) = a{p) + e{p) (mod p) 

Thus, if ere = cr(^?) x e{p) = —1 , then p | t/p+i = Up-cre ■ 

• With: m =p — 1 : 

2P-2[/p_^ = (p _ i)p?'-2 + _ i)£)pP-4 + ... + _ l)r>^p + 

2P-^Up-i{P^ -D) = -(p2)^ + i:)^ = £(p) _ (mod p) 

Thus, if ere = cr(p) X £{p) = 1 , then p \ C/p_i = Up-cre ■ 

□ 

Property (8.4.4) page 197 

1 — (T£ 

If p is an odd prime and p\ Q, then: Vp-ae = 2crQ 2 (mod p) . 
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Theorem 2 (8.4.1) If p is an odd prime and p \ QRD , then: 



( p\V 



p—cre 



when a = — r 



< 



p I U 



p—ae 
~2~ 



when a = T 



Definition (8.4.2) page 197 of a;(m) : For a given m, denote by a; = 
u){m) the value of the least positive integer k such that m \Uk . If u){m) 
exists, u}{rn) is called the rank of apparition of m . 

Theorem 3 (8.4.3) 



Theorem 4 (8.4.5) // (m, Q) = 1 , then u{m) exists. 

Theorem 5 (8.4.6) // (iV, 2QRD) = 1 andN±l is the rank of appari- 
tion of N, then N is a prime. 

Theorem 6 (8.4.7) // {N,2QRD) = 1 , C7iv±i = (mod AT) and 
U N±i 7^ (mod N) for each distinct prime divisor q of N ± 1, then N 



Proof: 

Let CO = u;{N) . We see that uj \ N ± 1 , hut uj \ (N ± l)/q . Thus if 
II iV ± 1 , then q" \ lu . It follows that w = ± 1 and A'^ is a prime by 
Theorem 5 (8.4.6) . 



3 Fn prime ^> F„ | V Fn-i and F„ | 52n_2 



If k \ n, then U ^ \ Un ■ 
If m I Un, then u{m) \ n . 



1 

is a prime. 
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Let N = Efi = 2^" + 1 with n > 1 be an odd prime. 
Let: P = ^/R,R = 7,Q = l,andD = P^-4Q = 3. 

Hereafter we compute ("^/n) and (Vn) • 

• (%r) • f ^ prime 



Since: < 



iV = (4)2" +1 = 2 (mod 3) 




= -1 . 




[ m = m X (-i)"^"^ 
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• (7/n) : We have: J 2^ =1 (mod 7) 

I 23"+'' = 2^ (mod 7) 

With 2" = 6 (mod 3) , we have: 2^" + 1 = 2'' + 1 (mod 7) . Then we study 
the exponents of 2, modulo 3 . We have: 2^ = 1 (mod 3) , and: 

22"» = 1 (mod 3) 
If n = 2m { N = 22^"" + 1 = 2^ + 1 = 3 (mod 7) 

m = m = -1 



If n = 2m + 1 < 



22m+i = 2 (mod 3) 
N = 22'"^+' + 1 = 22 + 1 = 5 (mod 7) 

I m = m = -1 



7-1 02*^ 



Finally, we have: (7/n) = (^/7)(-l)^^ = (%) 



(%) 


= m = 


-1 


(%) 


= m = 


-1 


(Q/n) 


= (Vn) = 


+1 



So we have: 



Since a = — r , ere = +1 , and F„ | Qi?!) with n > 1, then by Theorem 2 
(8.4.1) we have: 

Fn prime ^ F„ | V Fr,-i = K2"-i 



By (1) we have: V^k-i = 'S'fe_2 and thus, with k = 2"': F„ | S'2n_2 • 



□ 



4 Fn I S'2n_2 ^> F„ is a prime 

Let N = Fn with n > 1 . By (1) wc have: N \ S'2«-2 =^ ^ I ^2"-i ■ 
And thus, by (4.2.6) page 74 ( U2a = UaVa ) , we have: | U^2'^ . 

By (4.3.6) page 85: ( {Vn, Un) \ 2Q" for any n ), and since Q = 1 , then: 
(V^2"-i) ^22"-i) ~ ^ ^^"^ thus: N I L'"22"-i since N odd. 



With w = u{N) , by Theorem 3 (8.4.3) we have : a; | 2^" and lo f 22"-i . 

This implies: a; = 2^" = AT - 1 . Then A" - 1 
and thus by Theorem 5 (8.4.6) N is a prime. 



This implies: u = 2'^^ = N — 1 . Then AT — 1 is the rank of apparition of N, 

□ 
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This test of primality for Fermat numbers has been communicated to the 
community of number theorists working on this area on mersenneforum.org 
fhttp:/ /www.mersenneforum.org/showthread.php?t=2130 ) in May 2004, and 
the proof was finahzed in September 2004. 

Then, in a private communication, Robert Gerbicz provided a proof of the 
same theorem based on (5[-v/2T]. 



5 Numerical Examples 

(mod F2) 5*0 = 5 i 6 'S'22-2 = 

(mod F3) = 5 i 23 3> 13 3> 167 ^^ 131 ^^ 197 = -60 ^ Sgs^a = 
(mod F4) = 5 23 3> 527 3> 15579 ^^ 21728 42971 ^^ 1864 
1033 18495 3> 27420 ^ 15934 2016 ^ 960 34 4080 ^ 82^-2 = 



6 Appendix: Table of Ui and Vi 

With n = 2, 3, 4, we have the following (not proven) properties (modulo -F„,): 



Upn-b = 


5 


= 


-23 


Upn-i = 


6 


^F„-4 = 


-4 


Upn-i = 


1 


VPn-S = 


-5 


UFn-2 = 


1 


yFn-2 = 


-1 


Upn-l = 





ypn-i = 


-2 


Upn = 


-1 




-1 


U Fn+l = 


-1 


^Fn+l = 


-5 


U Fn+2 = 


-6 


VFn+2 = 


-4 


, UFn+3 = 


-5 


, ypn+S = 


-23 



The values of {/'„ and V'n (n>i) with (P, Q) = (\/3, —1) can be built by: 

U'2n = U2n \ ^'2n = ^2n 

f^'2n+l = V2n+l [ ^'2n+l = f^2n+l 

Values of Ui and l/j in previous tables can be computed easily by the following 
PARI/gp programs: 

J72j+i : U0=1;U1=6; for(i=l,N, U0=5*U1-U0; U1=5*U0-U1; print (4*i+l , " 
",U0); print (4*i+l," ",U1)) 
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Table 1: P =8v/7 , Q = l 



i 


Ui (mod Fi) 


(mod Fi) 
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1 


1 


1 


2 


1 
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4 
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6 


4 
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4 


1 


8 
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Table 2: P = ^/7 , Q = 1 , Modulo Fi 



i 


Ui (mod Fo) 


Vi (mod Fo) 








2 


1 


1 


1 


2 


1 


5 


3 


6 


4 


4 


5 


6 


5 


12 


2 


6 


7 


8 


7 


3 


6 


8 


13 





9 


3 


11 


10 


7 


9 


11 


12 


15 


12 


5 


11 


13 


6 


-4 


14 


1 


-5 


15 


1 


-1 


16 





-2 


17 


-1 


-1 


18 


-1 


-5 


19 


-6 


-4 


20 


-5 


11 


21 


5 


15 


22 


10 


9 


23 


14 


11 


24 


4 






Table 3: P = ^7 , Q = l , Modulo F2 
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i 


Ui (mod F3) 


Vi (modFa) 
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1 


1 


1 


2 


1 


5 
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4 
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167 
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64 


38 
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1 


-1 


256 





-2 
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-1 


-1 


258 


-1 


-5 


259 


-6 


-4 


260 


-5 


-23 



Table 4: P = ^7 , Q = l , Modulo F3 
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i 


Ui (mod F4) 


Vi (modF4) 


2048 


9933 


15934 


4096 


567 


2016 


8192 


28943 


960 


16384 


63129 


4080 


32768 


5910 





65532 


5 


-23 


65533 


6 


-4 


65534 


1 


-5 


65535 


1 


-1 


65536 





-2 


65537 


-1 


-1 


65538 


-1 


-5 


65539 


-6 


-4 


65540 


-5 


-23 



Table 5: P = ^7 , Q = l , Modulo F4 
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